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Negative heat capacity in the critical region of nuclear fragmentation: an 
experimental evidence of the liquid-gas phase transition 
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An experimental indication of negative heat capacity in excited nuclear systems is inferred from 
the event by event study of energy fluctuations in Au quasi-projectile sources formed in Au + Au 
collisions at 35 A.MeV. The excited source configuration is reconstructed through a calorimetric 
analysis of its de-excitation products. Fragment partitions show signs of a critical behavior at about 
5 A.MeV excitation energy. In the same energy range the heat capacity shows a negative branch 
providing a direct evidence of a first order liquid gas phase transition. 



Phase transitions are the prototype of a complex sys- 
tem behavior which goes beyond the simple sum of indi- 
vidual properties Q. In macroscopic systems the ther- 
mostatistical potential presents non analytical behaviors 
which unambiguously marks a phase transition. Non 
analytical behaviors of infinite systems originate from 
anomalies of the thermostatistical potentials in finite sys- 
tems ly,^. Specifically in microcanonical finite systems, 
the entropy is known to present a convex intruder in 1-st 
order phase transitions associated to a negative heat ca- 
pacity between two poles. A 2-nd order phase transition 
is characterized by the merging of the two poles. 

The experimental study of phase transitions in finite 
systems has recently attracted a strong interest from var- 
ious communities. Bose condensates with a small number 
of particles [Q , melting of solid atomic clusters g , vapor- 
ization of atomic nuclei H are examples of attempts to 
study phase transitions in finite systems. The problem 
usually encountered with these small systems is how to 
control the equilibrium and how to extract the thermo- 
statistical variables from observable quantities in order 
to identify the possible phase transition. This is for in- 
stance the case in heavy ion reactions in which excited 
nuclear systems are formed. Comparingthe observed de- 
cay channels with statistical models [^,0 it seems that a 
certain degree of equilibration is reached pM but up to 
now it has not been possible to unambiguously identify 
the presence of the expected liquid-gas phase transition. 

It has recently been shown [|| that for a given total 
energy the average partial energy stored in a part of the 
system is a good microcanonical thermometer while the 
associated fluctuations can be used to construct the heat 
capacity. In the case of a phase transition anomalously 
large fluctuations are expected as a consequence of the 
divergence and of the possible negative branch of the heat 
capacity. Let us consider an equilibrated system which 
can be decomposed into two independent components so 



that the energy is simply the sum of the two partial en- 
ergies Et — El + £'2 and that the total level density 
Wt = exp{St) is the folding product of the two partial 
level densities Wi = exp[Si). 

An example of such a decomposition is given by the ki- 
netic and the potential energies in the absence of velocity 
dependent interactions. 

The probability distribution of the partial energy 
Pf' [El) = Wi{Ei)W2{Et - Ei)/Wt{Et) is peaked on 
a most probable energy partition which requires the 
equality of the two microcanonical temperatures T~ = 
dSJdEi, i = 1, 2 of the two subsystems (Ti = T2 = T). 
The fluctuations |l| around a maximum Ei of P^ * are 
given by 
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where: C"'^ = -T^d'^Si/dEf {i = 1, 2) are the heat ca- 
pacities calculated for the most probable energy partition 
{El ,Et — El). 

The heat capacity C^^ = -T^d^St/dEf reads: 



Ct ~ Ci + C2 



CI 



Ci - al/T^ 



(2) 



In the absence of a phase transition because of the mi- 
crocanonical constraint one expects the fluctuations al 
to be smaller than the canonical expectation CiT"^ and 
the total heat capacity to be positive. Phase transitions 
are signed by poles and negative heat capacities corre- 
sponding to anomalously large fluctuations a\ > CiT^ . 

A word of caution is here necessary. Convex intrud- 
ers are expected for phase transitions of systems charac- 
terized by a unique extensive observable e.g. the total 
energy, as for example in the case of the melting transi- 
tion which is characterized by a symmetry breaking. In 
the liquid-gas case, since the volume is directly related to 
the order parameter it should be considered as another 



extensive variable. In the experimental situation, since 
the break up volume is at best known only in average, 
the pressure, interpreted as the Lagrange multiplier as- 
sociated with the volume observable, appears to be the 
relevant state variable together with the total deposited 
energy. In such a case the energy fluctuations are related 
to Cp and a convex intruder is expected. In principle 
one could argue that the same information on the heat 
capacity can be obtained by taking the derivative of the 
correlation between the temperature and the excitation 
energy (the so called caloric curve). We want to stress 
that fluctuations are a characteristic of the state and so 
depend on the pertinent state variable while the caloric 
curve T{E) depends upon the specific thermodynami- 
cal transformation from one state to another. Therefore, 
the information obtained by taking the derivative of the 
measured caloric curve may differ from the information 
coming from the fluctuations. 

Here we report on a study on the de-excitation proper- 
ties of Au quasi-projectiles formed within an excitation 
energy range from 1 to 8 A.MeV. The experiment was 
performed at the K1200-NSCL Cyclotron of the Michi- 
gan State University. Beams of Au ions at 35 A.MeV in- 
cident energy were used to bombard Au foils. The MUL- 
TICS and MINIBALL arrays were coupled to measure 
light charged particles and fragments with a geometric 
acceptance greater than 87% of 47r. For experimental 
details see ref. Q. Peripheral collisions of a predom- 
inantly binary character have been selected by requir- 
ing the velocity of the largest fragment in each event to 
be at least 75% of the beam velocity. For each event, 
the fragments were considered as originating from the 
quasi-projectile if forward emitted in the centre of mass 
reference frame. The contribution of light particles is 
then added. To avoid pollution from other sources than 
the quasi-projectile, the backward emitted particles are 
substituted by the symmetric of the forward emission in 
the quasi-projectile reference frame. Only events where 
the total source charge is within 10% of the original Au 
charge have been taken into account, to deal with com- 
parable finite size effects over the whole excitation energy 
range. The total excitation energy of the source is mea- 
sured from calorimetry on an event by event basis. This 
allows a sorting of the events as a function of energy as 
in the microcanonical statistical ensemble. 

A necessary condition for the application of statisti- 
cal equilibrium concepts is spherical symmetry in mo- 
mentum space, which can be verified by looking at the 
fragment angular distribution in the reaction plane. As 
an example Fig.l displays the fragment {Z > 3) angular 
distribution in the source reference frame for different in- 
tervals of the excitation energy. The observed isotropy 
in the source rest frame (only slightly affected in the 
backward hemisphere by the experimental acceptance) 
is a guarantee of the quality of the source reconstruction 
analysis |^. 

Indications of the possible occurrence of a phase tran- 
sition can be seen in the critical behavior of the fragment 



partitions. Specifically, in the vicinity of a second order 
phase transition (in infinite systems), partitions follow a 
finite size scaling characterized by two independent crit- 
ical exponents. This generic behavior allows to classify 
critical systems in various universality classes which de- 
pend upon few characteristics such as the dimensionality 
of the space and of the order parameter. In the same 
class, various systems present a similar phenomenology 
and in particular identical critical exponents. 

• EyA„=5-8 O r/A„=3-5* r/A„=1-3 (A.MeV) 
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FIG. 1. Angular distribution of the fragments (Z > 3) in 
the guasi-projectile reference frame in different intervals of the 
excitation energy per nucleon of the source E* /Aq. 

Experimentally, information on critical exponents can 
be inferred from the analysis of the Campi scatter plot 
which gives the correlation between the heaviest frag- 
ment produced in each event and the corresponding sec- 
ond moment of the charge distribution [|lO|,Q. This plot, 
presented in Fig. 2a) for our data, shows two branches 
characteristic of a subcritical regime and a supercriti- 
cal one. The zone where the two branches join corre- 
sponds to the critical region. The critical exponents can 
be extracted from the mean values of the moments of the 
Campi plot, giving P/j = 0.29±0.01. Fig. 2 b) || shows 
the charge distribution in the critical region of the Campi 
plot (events falling into the selected area). One can rec- 
ognize the critical distribution in the form of a power law 
of exponent r = 2.13 ± .08 in agreement within the errors 
with the value of /3/7, according to the scaling relation 
T — 2 — l3/"f/{l+(3/"f). Similar results were already found 
by the EOS collaboration (l|]. 

A qualitative indication that the observed phase tran- 
sition has a thermodynamical origin comes from the 
fact that, in the subcritical regime, an exponent /3 = 
0.33±0.04 can be obtained by fitting with a power law the 
charge of the heaviest fragment as a function of the exci- 
tation energy per nucleon E* / Ao of the quasi-projectile 
(Fig. 2c)), which means that E* /A{) can be considered as 
a critical parameter. 
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action, asymptotic kinetic energies have to be corrected 
from the Coulomb boost. To take into account these dis- 
tortions, primary partitions have been reconstructed by 
applying the energy balance at the production time event 
by event: 



mo + £■* = .Bi + .^2 



(3) 



where toq is the mass of the source, E* the excitation 
energy calculated via calorimetry, and the energy E2 can 
be chosen as: 



Eo 



M 



Er 



(4) 
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Here rrii is the mass excess of primary fragment i, M is 
the event multiplicity at the production time, and Ecoui 
is the Coulomb energy of the partition. 



FIG. 2. (color) a) Campi scatter plot: Logarithm of the 
charge of the largest fragment in each event (Zug) as a func- 
tion of the normalized second moment of the charge distri- 
bution, c) Logarithm of Zttg as a function of the excitation 
energy of the source, h) Charge distribution for the critical 
region of the Campi scatter plot (dashed contours of panels a) 
and c)). 

The solid symbols represent the mean correlation. The full 
lines are fits resulting from power-law behaviors. 

The extracted values of the critical exponents are per- 
fectly compatible with a liquid-gas transition but also 
with a geometrical percolation type of transition. More- 
over, as found in the Lattice Gas Model pi] for finite 
systems, a critical behavior in fragment observables can 
also be consistent with a phase coexistence of a 1-st or- 
der phase transition. Indeed in the coexistence region of 
small systems the liquid cluster is not much larger than 
the vapor fragments and may mimic critical fluctuations 
of the mass distribution with approximately the same 
critical exponents as at the critical point. In other words 
the analysis of the charge distribution in terms of critical 
behavior and critical exponents is a sign of a phase tran- 
sition but does not allow to determine the corresponding 
universality class nor the associated order. Therefore, a 
more direct information about the order of the transition 
is mandatory. 

To further progress on this point we shall now examine 
the partial energy fluctuations which are a direct observ- 
able to explicitly test thermodynamical equilibrium and 
extract thermodynamical state variables from the exper- 
imental data. In order to exploit eq.(|2|), one has to find a 
suitable decomposition of the total energy Et into Ei and 
E2 ■ In the case of nuclear fragmentation data, the repar- 
tition of the total energy is complicated by the fact that 
thermal equilibrium fragments are produced hot, while 
the detected fragments are collected at infinity after sec- 
ondary de-excitation, i.e. with lower mass. Moreover, 
because of the presence of the long range Coulomb inter- 
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FIG. 3. (color) a),c) Partial energy per nucleon Ei/Ao as 
a function of the excitation energy E* /Aq. The full line gives 
the average values (Ei/Ao). b), d) Normalized variance of 
Ei/Ao as a function of E* /Aq. Left (right) panels refer to 
the freeze-out hypotheses I (II) (see text). 

The primary masses have been obtained by sharing 
the charge detected in the form of light particles and 
neutrons among the detected fragments, following two 
extreme freeze-out hypotheses [^,0 . In the first assump- 
tion (/) primary fragments {Z > 3) have a charge-to- 
mass ratio as in the entrance channel, so that they de- 
excite only through neutron evaporation and all the light 
charged particles are considered as primary. In the sec- 
ond hypothesis (//) the totality of both neutrons and 
light charged particles is shared among final fragments. 

The experimental correlation between the partial ener- 
gy Ei/Aq and the excitation energy per nucleon of the 
source E* /Aq is shown in Fig. 3a), c) for the two freeze- 
out hypotheses. The normalized variance of Ei/Aq, in 
bins of E* /Aq, is represented in Fig. 3b), d). The two 



freeze-out hypotheses are in quahtative agreement, and 
the presence of peaks of Fig. 3 b, d) indicate a phase tran- 
sition governed by the equihbrium between the kinetic 
and potential degrees of freedom ||3|], i.e. a thermodynam- 
ical phase transition, in both cases. FinaUy, the actual 
microcanonical temperature T can be obtained from (Ei) 
by inverting the kinetic equation of state 



M 



{E^)^{J2a.)T' + {^{M-l))T 



(5) 



i=l 



where Ui is the fragment dependent level density param- 
eter. The brackets < > indicate the average on the 
events with the same E* . This also allows the extraction 
of the kinetic heat capacity Ci obtained by taking the 
numerical derivative of (Ei) with respect to T. 



It is important to remark that the thermodynami- 
cal transition energy interval corresponds to the criti- 
cal region of the Campi scatter plot: the observation of 
power laws and critical behaviors in the fragment charge 
distribution is directly correlated to the occurrence of 
anomalous fluctuations indicating a negative heat capac- 
ity. These two pieces of informations which correspond 
to the theoretical expectations |lj,^,|l5[ indicate that a 
thermodynamical first order phase transition has taken 
place in finite nuclear systems formed in nuclear reac- 
tions. 
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FIG. 4. Specific heat capacity (solid symbols) obtained from 
eq.(2). The panel on the left (right) refers to the freeze-out hy- 
potheses I (II). The grey contour indicates the Ct distribution 
from the combination of the different uncertainties. 

Fig. 4 shows Ct from eq.(0). The Ct distribution in 
Fig. 4 takes into account the uncertainties in the numeri- 
cal derivatives of eq.(||) as well as statistical errors. 

A clear negative branch appears pointing to a 1-st or- 
der liquid-gas phase transition, the distance between the 
poles being associated with the latent heat. From Fig.s 3 
and 4 it is also apparent that the freeze-out hypothe- 
sis is not crucial for our main conclusions, if we do not 
pretend to give a quantitative evaluation of the latent 
heat. The persistency of the signal has been checked 
against variations of the excitation energy bin and the 
different hypotheses on the energy reconstructions such 
as the number and energy of the (undetected) neutrons. 

This phenomenon predicted as a specific feature of 1- 
st order phase transitions in finite systems in the mi- 
crocanonical ensemble |lq] to our knowledge was never 
experimentally observed insofar. 
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